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Abstract 

The production of a Biggs boson in association with a VP-boson is the most 
likely process for the discovery of a light Biggs at the Fermilab Tevatron. Since 
it decays primarily to 6-quark pairs, the principal background for this associated 
Biggs production process is Wbb, where the bb pair comes from the splitting of an off 
mass shell gluon. In this paper we investigate whether the spin angular correlations 
of the final state particles can be used to separate the Biggs signal from the Wbb 
background. We develop a general numerical technique which allows one to find a 
spin basis optimized according to a given criterion, and also give a new algorithm 
for reconstructing the W longitudinal momentum which is suitable for the WH and 
Wbb processes. 



1 Introduction 



At present, the existence of a neutral Higgs boson is certainly the largest unresolved 
problem in the standard model (SM). Its mass is a priori unknown, but direct searches 
and precision electroweak measurements constrain it to be 90 < Mu < 280 GeV at 95% 
confidence level M. At the Tevatron collider there is a possibility to search for the SM 
Higgs using the decay mode H — > bb @, and the most promising process is the associated 
Higgs production 

pp -> W(-> eu)H(^ bb) . (1) 

The Fermilab search is extremely important, especially because the mass range which can 
be covered at the Tevatron (100 < M# < 130 GeV) is also one of the most challenging 
regions for the LHC to look for the SM Higgs p[. With sufficiently large data sample the 
Higgs signal could be extracted from the background by analyzing the bb mass distribution. 
However, given the fact that there are several large backgrounds to process ([!]), any 
technique which can provide additional handles on distinguishing the signal from the 
background would be useful. 

In this paper we investigate the possibility of using the spin angular correlations for 
separating the associated Higgs production from its principal background at the Tevatron, 
the Wbb process 

pp-> W(-* ev)g*{-^ bb) . (2) 

In the case of e + e~ — *> ZH/ZZ in Q it was shown that spin angular correlations can pro- 
vide useful information if good spin bases are chosen. Since the qq — > WH/Wbb processes 
have the same spin structure, it is natural for one to ask a question whether a similar 
analysis would be useful for distinguishing ([[]) and (0) at the Tevatron. However, due to 
the hadronic collider environment, and also to the complexity of the Wbb amplitudes, it 
is obvious that in this case a numerical approach for finding the best spin basis is more 
appropriate than the approach used in |4| . For that reason we develop here a new method 
which allows one to find a spin basis optimized according to a given criterion. This tech- 
nique is completely general in the sense that it can be used for optimizing spin basis 
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regardless of which or how many processes are being considered. We apply our method 
to WH and Wbb processes, and suggest several possible strategies which could add new 
information in an experimental analysis. We also discuss one of the major uncertainties 
related to our analysis, and that is the W momentum reconstruction. Our results indicate 
that the method which has been used in the literature can distort angular distributions 
considerably, and is therefore inadequate for our purposes. Because of that we propose a 
new W reconstruction algorithm whose effects on angular distributions are significantly 
less destructive. 

The remainder of the paper is organized as follows: in Section ^| we give all relevant 
definitions, describe numerical method and suggest possible strategies for finding the 
optimal spin basis. In Section |3| we present our results for angular distributions, and show 
the effects which W reconstruction algorithm has on those. Conclusions are contained in 
Section |j. 

2 Angular correlations 

In order to apply the generalized spin-basis analysis || to processes ([I]) and (§) we first 
define the zero momentum frame (ZMF) production angle 6* (0 < 6* < it) as the angle 
between the incoming up-quark and the VT-boson produced in the qq' — > WX process (see 
Figure [IJ), where X is either H or g*. The spin states for W are defined in its rest frame, 
where we decompose its spin along the vector sw, which makes an angle £ with the X 
particle momentum in the clockwise direction. The X particle's spin can be decomposed 
in a similar way. Relationship between £ and 9* determines specific spin basis in which 
one can calculate angular correlations among the WH and Wg* decay products in ([]]) and 
(H). These correlations involve distributions of the angle \w (Xbl) that the charged lepton 
(6-quark) makes with the spin vector of the VT-boson (bb system) . Figure illustrates the 
definitions for angles £ and xw- 

In the case of e + e~ — ► ZH/ZZ — > // jets the procedure for finding the optimal 
spin basis was based on separating the polarized amplitudes for e + e~ — * ZH/ZZ M. In 
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particular, it was shown that very good separation between the ZH and ZZ events can be 
obtained in the transverse basis, in which the longitudinal component of the ZH matrix 
element is zero by construction. 

Since the amplitude for the process qq' — > WH has the same spin structure as the 
one for e + e~ — > ZH,^ the transverse basis is also a good starting point for examining the 
cos x distributions in the WH and Wbb processes. It is defined by 

tan#* 

tan £ = / > ( 3 ) 



1-^ 

where f3\y is the ZMF speed of the H^-boson. Nevertheless, due to the complex nature of 
the qq' — > Wg* amplitudes, and also to the fact that in pp collisions the center-of-mass 
energy Vs is not fixed, the approach of Ref. M for finding the optimal spin basis is not 
practical for our purposes here. Because of that, instead of trying to separate polarized 
cross sections for qq' — > WH/Wg*, we attempt to find the best basis for processes (]l|) and 
(0) by distinguishing the cosx distributions directly, using a suitable multidimensional 
maximization procedure. 

The basic idea of our method is to divide cos 9* — cos £ plane into n x m regions, and 
to associate with each of those a histogram containing distribution in cosx-0 A specific 
spin basis is defined by choosing one of the cos£ bins for all of n cos 9* bins, while the 
total cos x distribution is obtained by summing contributions over the entire cos 9* range. 
In other words, if cos £j describes the spin vector of W (or bb system) in the i-th cos 9* 
bin, and <jj is the corresponding contribution to the cross section, we have 

da _ dai(cos£i) 
dcc&x dcosx 

In this way, by changing the n cos£j variables using multidimensional maximization al- 
gorithm, one can easily vary the definition of the spin basis until the optimal separation 

1 The spin structure for the process ud — ► WH can be found in Eqs. (4)-(6) in [Q. The full matrix 
element squared, including the decay of W^-boson, can be obtained from Eq. (2) in 

2 In general we do not know the up-quark momentum direction. However, the up-quark comes from 

the proton beam more than 95% of the time at the Tevatron, and therefore we will use proton direction 
instead of the up-quark direction in defining cos 9* for the rest of this paper. 
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of the WH and Wbb events is achieved. 

Results of this procedure will clearly depend on which criterion is used for determining 
the best possible separation of the signal and the background. We investigate here two 
possible criteria. The first one is based on distinguishing between the shapes of the cos x 
distributions for the two processes, and the function which we decided to maximize is 
given by 



f , da WH 1 r 

/ ctcosx^ cosx — / acosx 

a WH J a cos x (Jvtns. J 



d<J Wbh 

cosx 



a Wb iJ dcosx 



(5) 



With this criterion the resulting spin basis tends to give cosx distributions which are 
asymmetric for the WH signal events and symmetric for the Wbb background events. 

The second criterion which we examine is based on maximizing significance S/ v~B, 
where S and B correspond to the number of events for the signal and background, re- 
spectively: 

rcos Xmax da WH 

o cc dcosx-; ? (o) 

J COS Xrnin « COS \ 

rcos Xma da 

B oc J dcosx~. • v) 

J COS Xrnin U, COS X 

Once a particular spin basis is chosen and the cosx distribution for both processes is 
calculated using (|j), we choose angles Xmin and Xmax in such a way to maximize the ratio 
S/ \J~B. Note that in the above coefficients of proportionality include the NLO i^-factors, 
our assumptions on the integrated luminosity and double 6-tagging efficiency, etc. 

The main advantage of the method described above is that it offers a systematic ap- 
proach for investigating the possibility of using the spin angular correlations to distinguish 
signal events from the background, regardless of which or how many processes are being 
considered. For example, even though we are concerned here only with the leading or- 
der Wbb process as the most important background for the associated Higgs production 
at the Tevatron, it would be straightforward to include other backgrounds or the next- 
to-leading order effects as well. Note however that calculation of angular correlations 
between the spin vector of an intermediate gauge boson and momenta of its decay prod- 
ucts requires complete reconstruction of an event. That is a major difficulty in the case 
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of the WH/Wbb production where longitudinal component of neutrino is unknown. This 
issue will be discussed in more details in the following section. 



3 Numerical results 

Since the procedure outlined above requires large statistics in order to make errors in 
da jd cos x distributions as small as possible, for the results presented in this paper we 
generated about 10 8 events (for each process), using the VEGAS algorithm |7|]Fj Calcu- 
lations were done with n = 10 bins along the cos#* axis, and m = 1000 bins along the 
cos £ axis, while the search for the optimal basis was performed using the downhill simplex 



method [|H], 0. 

Even though the analysis described in the previous section can be performed for both 
W and bb sides of an event, we focus here only on the cosxw distributions. The reason 
is that the correlations on the W side of an event are much stronger and provide us with 
more distinguishing power for separating the WH and Wbb processes. [] 

All results shown in this paper are obtained for the W + production in pp collisions 
at VS = 2 TeV, with the MRSR1 parton distribution functions (a s (M z ) = 0.113) [0]. 
In order to improve our lowest order cross sections, instead of natural scales (/i m Mh) 



we used somewhat lower scale of ^ = 50 GeV [13|]. At this scale the NLO i^-factors are 
about 1.1 for both WH and Wbb processes. The Higgs mass was set to Mh = 120 GeV 
and the corresponding bb mass range to 102 < M b i < 141 GeV. In addition, we applied 



3 Because of the large statistics and the large number of histograms required by our method, Monte 

Carlo simulations which would include all other background processes, or take into account next-to-leading 

order corrections, would have to be done using a parallel event generator ||, JjJ. 

4 On the bb side of the event we were unable to find a spin basis which would considerably improve the 

small difference between the WH and Wbb processes that was obtained using the helicity basis. 
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the following set of isolation cuts and cuts on rapidity and transverse momentum: 



Rbbi Reb, R e b 


> 


0.7, 


\Vb\, \m\ 


< 


2, 


\Ve\ 


< 


2.5, 


\Pbl \pI\ 


> 


15 GeV 


\Pel \PI\ 


> 


20 GeV 



(8) 



Note the \p^\ cut is the missing Et cut and that the above cuts do not include a cut 
on cos 9* [ [n{] . Our results indicate that imposing the cos^* cut actually worsens our 
ability to separate the two processes based on the shape of their cos xw distributions, and 
therefore we did not include it in the simulations based on maximization of Eq. (|5]). On 
the other hand, it is well known that this cut can improve the S/ y/~B ratio by about 10% 



njj. Because of that, we take it into account for simulations based on the significance 
criterion. 

We first discuss our results obtained with the shape criterion. In this case we found 
that the optimal basis can be well approximated by the polynomial of the form 

k 

cos^^a^cosfl*) 2 ^ 1 . (9) 

i=i 

In particular, in Figure^ we compare the exact result obtained by maximization procedure 
to polynomial with k = 3 and coefficients 

ai = 0.2354 , 

a 2 = 0.1808 , (10) 
a 3 = -1.442 . 

There we also show the transverse basis with specific choice of f3y/ = 0.67, which is 
close to the averages of the distributions for both WH and Wbb processes (0.68 
and 0.66, respectively). The actual normalized cosxw distributions corresponding to the 
polynomial approximation of the optimal basis and for the transverse basis are given in 
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Figure [|. As expected, in the optimal basis the Wbb distribution is nearly symmetric. 
Figures ^ and |^ illustrate what one might expect in terms of the number of events per 
bin in those two bases. These results were obtained by multiplying our W + cross sections 
by four to take into account contributions from the W~ production and the contribution 
from the W ± decays into muons, by taking into account the NLO i^-factor of 1.1 for both 
WH and Wbb processes, and also by assuming the double 6-tagging efficiency of e% = 0.45 
and integrated luminosity of 10 fb _1 . We would like to point out here that the shape of 
the Wbb cos xw distribution is significantly different in the two bases being discussed. On 
the other hand, this is not the case for the WH process. Clearly, the difference in the 
shape of the cos xw distributions under the change of spin basis may provide an additional 
handle for separating the two processes. 

Another interesting possibility of using angular correlations for distinguishing between 
the signal and the background is illustrated in Figures [7] through |9|. Instead of looking at 
da/dcosxw directly, we investigate the M b i distributions of the quantity crcosxw- Those 
distributions vanish in the spin basis in which da jd cos xw is perfectly symmetric, because 
in evaluating d(a cos xw) I dM h i one effectively integrates over cos xw- This is precisely the 
reason why our optimal basis reduces the background in acosxw much more efficiently 
than the transverse basis, as can be seen by comparing Figures || and However, at this 
point one should also observe that the main disadvantage of analyzing quantities such as 
a cos xw is the inclusion of the statistical errors of the entire cos xw distribution, which 
may limit its potential usefulness in an experimental analysis with small statistics. 

In our simulations based on maximizing significance the number of events for both 
signal and background was obtained by summing all decay channels, and under the same 
assumptions as before (X-factors of 1.1, = 0.45, and / C dt = 10 fb _1 ). As mentioned 
earlier, besides cuts given in (|S|), here we also take into account a cut on cos 9* JTjJ], 



| cos 0*| <cos9* max , (11) 

which can increase the ratio S/\fB by about 10% (see Table [l|). Using the method 
described in Section ^|, and for any given value of cos^^, we were able to find a spin 



S 



basis (and a set of cuts on cos x) in which one could further improve this ratio by additional 
2-3%. 

Our results with cos0 max = 0.8 are shown in Figures [lCj through [13|. Figure [U] shows 



the optimal basis definition. In this case we found that it can be approximated by 

cosf = -0.857 sign(cos0*) + 0.391 cos 0* . (12) 

Normalized cos xw distributions corresponding to the optimal basis are compared in Fig- 
ure |ll| to the results obtained using the transverse basis. Figures [12] and |TJ illustrate 
what can be expected in terms of the number of events per bin in those two bases. 

Due to the fact that the longitudinal momentum of the neutrino is unknown, recon- 
struction of an event involving W^-boson is the most important problem related to the 
calculation of the spin angular correlations which we discussed in this paper. By assuming 
that W is on shell, and using p e and which are actually measured, this component can 
be reconstructed up to a two-fold ambiguity for a solution of a quadratic equation. The 



algorithm for choosing the correct solution which has been used in the literature [FJJ] is 



based on the asymmetry of the neutrino rapidity distribution. From Figure |T4] it can be 



readily seen that by choosing the larger (smaller) solution for p z y in the case of W + (W~), 
one can improve the probability of finding the correct W momentum. Nevertheless, we 
propose here that for the WH and Wbb processes reconstruction algorithm is based on 
the distribution of the difference between the W rapidity and the rapidity of the bb sys- 



tem (see Figure [15]). Since this distribution is peaked at zero, our prescription consists of 
choosing the solution for p z v which results in a smaller absolute value for r]w — 77^. The 
advantage of using rjw — Vbb instead of r] u is that its distribution is narrower. Furthermore, 
unlike the r] u distribution, it is almost identical for WH and Wbb, which means that our 
algorithm will work equally well for both processes. 

In order to investigate the effect that W reconstruction algorithm has on cosxw 
distributions, we have repeated calculations shown on Figure |] (without cuts on cos0*) 
for the polynomial approximation of the optimal basis (Eqs. (§) and (]T0|)), and for the 
transverse basis. Results given in Figures [TO through [19] show that the cos xw distributions 
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obtained using our prescription for reconstructing the W momentum are much more closer 
to the exact curves than are the ones obtained using the r] v algorithm. Note that one 
of the reasons for distortion of the reconstructed cosxw distributions is the fact that in 
our calculations the W width is taken into account, while the reconstruction algorithms 
assume the on-shell W. 

Besides the issues related to reconstruction of the W momentum, another problem 
which might affect experimental analysis of the cosxw distributions is the mismeasure- 
ment of the 6-quark momenta. We have simulated that by imposing a Gaussian distribu- 
tion of relative errors (with the variance of 5%) on both b and b momenta, and our results 
indicate that these effects are small. 

4 Conclusions 

In this paper we investigated the possibility of using the spin angular correlations for 
distinguishing between the WH and Wbb processes at the Fermilab Tevatron. We devel- 
oped a general numerical method for finding the spin basis optimized according to a given 
criterion, and also suggested several possible strategies for utilizing this technique in the 
Higgs search at the Tevatron. 

Our simulations indicate that the spin angular correlations may provide additional 
handle on separating the signal from the background. Still, there are several problems 
that would have to be solved for a successful experimental analysis, and the largest one 
is certainly the event reconstruction. In this regard we proposed a new W reconstruction 
algorithm which significantly reduces effects related to the W momentum ambiguities. 
We hope that this algorithm can be further improved upon. 

The obvious extension of this work would involve including the NLO corrections, as 
well including the other background processes. However, these calculations would be 
numerically quite challenging, and before they are attempted a feasibility study of their 
usefulness should be completed. 
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Table 1: Expected number of events at 10 fb for the signal (WH) and the background 
(Wbb) as a function of the cos 9* cut. Results shown are obtained for a 120 GeV Higgs. 



COS #max 


S 


B 


S/y/B 


1.0 


75 


260 


4.65 


0.9 


70 


198 


4.97 


0.8 


65 


161 


5.12 
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Figure 1: Scattering angle 9* in the zero momentum frame. The up-quark comes from 
the proton beam more than 95% of the time at the Fermilab Tevatron. 
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Figure 3: The optimal basis for the shape criterion (points), together with its k = 3 
polynomial approximation (solid line) and with the transverse basis for f3w = 0.67 (dashed 
line). 
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Figure 4: Normalized cos xw distributions for the polynomial approximation of the opti- 
mal basis (solid lines), and for the transverse basis (dashed lines). 



17 



50 



d 
•1-1 



40 — 



d 

> 



30 — 



Wbb 



a 

d 



20 — 



10 — 



WH 







J I I L 



J I I I I I I I L 



J I I L 



1.0 



-0.5 



0.0 
cosxw 



0.5 



1.0 



Figure 5: Distribution of the number of events per bin in the polynomial approximation 
of the optimal basis. The total number of events for WH is 75, and for Wbb is 260. 
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Figure 6: Distribution of the number of events per bin in the transverse basis. The total 
number of events for WH is 75, and for Wbb is 260. 
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Figure 7: M b i distribution of the cross section. No smearing of the 6-quark jet energies 
has been performed here. 
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Figure 8: M b i distribution of acosxw in the polynomial approximation of the optimal 
basis. 
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Figure 9: M b i distribution of a cos xw in the transverse basis. 
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Figure 10: The optimal basis for the significance criterion with cos6 max = 0.8 (points) 
together with its approximation given by Eq. fll2|) (solid line). 



23 



1.50 



1.25 — 



5 1-00 



Optimal basis 
Transverse basis 



/ - 




Figure 11: Normalized cosxw distributions for the basis optimized according to the 
significance criterion with cos9 max = 0.8 (solid lines), and corresponding results for the 
transverse basis (dashed lines). 
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Figure 12: Distribution of the number of events per bin in the optimal basis (cos max = 
0.8). Without cuts on cosxw we have 65 events for WH, and 161 events for Wbb 
(S/y/B = 5.12). With cuts —0.6 < cosxw < 1-0 these numbers are reduced to 60 
and 131, respectively {S/y/B = 5.24). 
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Figure 13: Distribution of the number of events per bin in the transverse basis {cos8 max = 
0.8). Without cuts on cosxw we have 65 events for WH, and 161 events for Wbb 
(S/y/B = 5.12). With cuts —0.1 < cosxw < 1-0 these numbers are reduced to 54 
and 110, respectively (S/y/B = 5.15). 
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Figure 16: Normalized WH cosxw distribution for the basis optimized according to the 
shape criterion. The true result is shown with the full lines, while results obtained using 
the rjw — Vbl an d f]v W reconstruction algorithms are plotted with the dashed lines and 
short dashed lines, respectively. 
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Figure 17: Normalized Wbb cosxw distribution for the basis optimized according to the 
shape criterion. The true result is shown with the full lines, while results obtained using 
the rjw — rj b i and r\ v W reconstruction algorithms are plotted with the dashed lines and 
short dashed lines, respectively. 
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Figure 18: Normalized WH cos xw distribution for the transverse basis. The true result is 
shown with the full lines, while results obtained using the rjw—Tjbb an d f} v W reconstruction 
algorithms are plotted with the dashed lines and short dashed lines, respectively. 
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Figure 19: Normalized Wbb cosxw distribution for the transverse basis. The true result is 
shown with the full lines, while results obtained using the r/w—Vbi an d t\ v W reconstruction 
algorithms are plotted with the dashed lines and short dashed lines, respectively. 
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